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We have studied the centrality dependence of charged particle forward-backward multiplicity 
correlation strength in Au+Au collisions at yijvjv=200 GeV with a parton and hadron cascade 
model, PACIAE, based on PYTHIA. The calculated results are compared with the STAR data. The 
experimentally observed correlation strength characters: (1) the approximately fiat pseudorapidity 
dependence in central collisions and (2) the monotonous decrease with decreasing centrality are 
well reproduced. However the theoretical results are larger than the STAR data for the peripheral 
collisions. A discussion is given for the comparison among the different models and STAR data. 
A prediction for the forward-backward multiplicity correlation in Pb+Pb collisions at y / ijviv=5500 
GeV is also given. 
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I. INTRODUCTION 

The study of fluctuations and correlations has been 
suggested as a useful means for revealing the mechanism 
of particle production and Quark-Gluon-Plasma (QGP) 
formation in relativistic heavy ion collisions P, [2J • Cor- 
relations and fluctuations of the thermodynamic quan- 
tities and/or the produced particle distributions may be 
significantly altered when the system undergoes phase 
transition from hadronic matter to quark-gluon matter 
because of the very different degrees of freedom between 
two matters. 

The experimental study of fluctuations and correla- 
tions becomes a hot topic in relativistic heavy ion colli- 
sions with the availability of high multiplicity event-by- 
event measurements at the CERN-SPS and BNL-RHIC. 
An abundant experimental data have been reported [3- 
Q where a lot of new physics arise and are urgent to 
be studied. A lot of theoretical investigations have been 
reported as well [7rjT3|. 

Recently STAR collaboration have measured the 
charged particle forward-backward (FB) multiplicity 
correlation strength 6 in a given centrality bin size at 
different centralities from the most central to peripheral 
Au+Au collisions at ^snn— 200 GeV ja|. The outstand- 
ing features of STAR data are: 

1. In most central collisions, the correlation strength 
is approximately independent of the distance be- 
tween the centers of forward and backward pseu- 
dorapidity bins, A77. 

2. The correlation strength monotonously decreases 
with decreasing centrality. 



3. In the peripheral collisions, the correlation 
strength approaches to an exponential function of 
A77. 

A lot of theoretical interests has been stimulated. 

The wounded nucleon model was used in [l2| to study 
the correlation strength, the first two characters of STAR 
data were reproduced, but the third one was not. In Ref. 

the Glauber Monte Carlo code (GMC) with a "toy" 
wounded-nucleon model and the Hadron-String Dynam- 
ics (HSD) transport approach have been used to analyze 
the STAR data. They used three different centrality de- 
terminations: the impact parameter hi, the number of 
participant (wounded) nucleons N part , and the charged 
particle multiplicity N c h in midrapidity \rj\ < 1. The 
first two characters of STAR data can be reproduced by 
the N part and N c h centrality determinations, while the 
third character can not. 

We have used a parton and hadron cascade model, PA- 
CIAE, to investigate the centrality bin size dependence 
of charged particle multiplicity correlation in most cen- 
tral (5, 0-5, and 0-10%) Au+Au collisions at ^5^=200 
GeV [lj]- It turned out that the correlation strength 
increases with increasing bin size. In fl4| we follow [15l ] 
defining the charged particle FB multiplicity correlation 
strength b as 



{n f n b ) 



(n f )(n b ) 



W) 2 



ov(nf, n b ) 
var(n / ) 



(1) 



'Corresponding author: sabh@ciae.ac.cn 



where n/ and n b are, respectively, the number of charged 
particles in forward and backward pseudorapidity bins 
defined relatively and symmetrically to a given pseudora- 
pidity r/. (n/) refers to the mean value of nt for instance. 
cov(rif,n b ) and var(n/) are the FB multiplicity covari- 
ance and forward multiplicity variance, respectively. 

In this paper we use the PACIAE model to study 
the charged particle FB multiplicity correlation strength 
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b in a given centrality bin size at different centralities 
from the most central to peripheral Au+Au collisions 
at v /s]v]v=200 GeV. A discussion is given for the com- 
parison among the models and STAR data as well as 
the STAR's convention of different centrality determina- 
tion for different measured At] points. A prediction for 
the forward-backward multiplicity correlation in Pb+Pb 
collisions at ■ s /snn=5500 GeV is also given. 



II. THE PACIAE MODEL 

The parton and hadron cascade model, PACIAE [l6| , 
is based on PYTHIA |l7[ which is a model for hadron- 
hadron (hh) collisions. The PACIAE model is composed 
of four stages: the parton initialization, parton evolu- 
tion (rescattering) , hadronization, and hadron evolution 
(rescattering) . 

1. PARTON INITIALIZATION: In the PACIAE 
model, a nucleus-nucleus collision is decomposed into 
the nucleon-nucleon collisions based on the collision ge- 
ometry. A nucleon-nucleon (NN) collision is described 
with the PYTHIA model, where a NN (hadron-hadron, 
hh) collision is decomposed into the parton-parton col- 
lisions. The hard parton-parton collision is described 
by the lowest-leading-order (LO) pQCD parton-parton 
cross section [l8[ with modification of parton distribu- 
tion function in the nucleon. And the soft parton-parton 
interaction is considered empirically. The semihard, be- 
tween hard and soft, QCD 2 — !• 2 processes are also 
involved in the PYTHIA (PACIAE) model. Because 
of the initial- and final-state QCD radiation added to 
the parton-parton collision processes, the PYTHIA (PA- 
CIAE) model generates a partonic multijet event for a 
NN (hh) collision. That is followed, in the PYTHIA 
model, by the string-based fragmentation scheme (Lund 
string model and/or Independent Fragmentation model), 
thus a hadronic final state is reached for a NN (hh) 
collision. However, in the PACIAE model the above 
fragmentation is switched off temporarily, so the result 
is a partonic multijet event (composed of quark pairs, 
diquark pairs and gluons) instead of a hadronic final 
state. If the diquarks (anti-diquarks) are split forcibly 
into quarks (anti-quarks) randomly, the consequence of 
a NN (hh) collision is its initial partonic state composed 
of quarks, anti-quarks, and gluons. 

2. PARTON EVOLUTION: The next stage in the 
PACIAE model is the parton evolution (parton rescat- 
tering). Here the 2 — s- 2 LO-pQCD differential cross 
sections [l8[ are employed. The differential cross section 
of a subprocess ij — > kl reads 
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where the factor K is introduced considering the higher 
order pQCD and non-perturbative QCD corrections as 
usual, a s stands for the strong (running) coupling con- 



stant. Taking the process q\q2 — > 9i?2 as an example 
one has 
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where s, t, and u are the Mandelstam variables. Since it 
diverges at t—0 one has to regularize it with the parton 
color screen mass /i 
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The total cross section of the parton collision, i + j, 
then reads 
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With the total and differential cross sections above the 
parton evolution (parton rescattering) can be simulated 
by the Monte Carlo method. 

3. HADRONIZATION: The parton evolution stage is 
followed by the hadronization at the moment of partonic 
freeze-out (no any more parton collision). In the PA- 
CIAE model, the Lund string fragmentation model and 
phcnomenological coalescence model are supplied for the 
hadronization of partons after rescattering. The Lund 
string fragmentation model is adopted in this paper. We 
refer to [16| for the details of the hadronization stage. 

4. HADRON EVOLUTION: After hadronization the 
rescattering among produced hadrons is dealt with the 
usual two-body collision model. Only the rescattering 
among 7r, k,p, n, p(u>), A, A, E, S, f2, J/\& and their an- 
tiparticles are considered in the calculations. An isospin 
averaged parametrization formula is used for the hh cross 
section ■ We also provide an option for the con- 
stant total, elastic, and inelastic cross sections (a 



NN 



N 



40 mb, o&t 



35 mb, ofJJ = 10 mb) and 



25 mb, a™ 

the assumed ratio of inelastic to total cross section of 
0.85. More details of hadronic rescattering can see [2l| . 



III. CALCULATIONS AND RESULTS 

In our calculations the default values given in the 
PYTHIA model are adopted for all model parameters 
except the parameters K and b s (in the Lund string frag- 
mentation function). K=3 is assumed and b s =6 is tuned 
to the PHOBOS data of charged particle multiplicity in 
0-6% most central Au+Au collisions at y / s]vjv=200 GeV 
[iH, as shown in Tab. |U 

In the theoretical calculation it is convenient to define 
the centrality by impact parameter bi. The mapping 
relation between centrality definition in theory and ex- 
periment 
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TABLE I: Total charged particle multiplicity in three r\ 
fiducial ranges in 0-6% most central Au+Au collisions at 
,/snn=200 GeV. 



AT c h(M < 4.7) N ch (\ V \ < 5.4) AT ch (total) 



PHOBOS" 
PACIAE 



4810 ± 240 
4819 



4960 ± 250 
4983 



5060 ± 250 
5100 



The experimental data are taken from [22]. 



is introduced [23J . In the above equation, g stands for the 
geometrical (total) cross section percentage (or charged 
multiplicity percentage) used in the experimental deter- 
mination of centrality. Ra = 1.12A 1 / 3 + 0.45 fm is the 
radius of nucleus A. We have also used the centrality 
determination based on the charged multiplicity N c h in 
the central pseudorapidity window \rj\ < 1 to study the 
FB multiplicity correlation strength as a function of A77. 
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FIG. 1: (Color online) The charged particle pseudorapidity 
distribution for the specified centrality bin in Au+Au colli- 
sions at t/Fnn=200 GeV. The solid and open symbols are 
the PHOBOS data [H] and PACIAE results, respectively. 

We compare the calculated charged particle pseudora- 
pidity distribution in specified centrality bin in Au+Au 
collisions at y/Jjv r /v=200 GeV with the corresponding 
PHOBOS data Q in Fig. 1. Here the theoretical cen- 
trality are defined by impact parameter bi via Eq.([B]). 
One sees that the PHOBOS data are well reproduced. 

The PACIAE model results of FB multiplicity corre- 
lation strength are shown in Fig. 2 (a) and (b) for the 
centrality determinations of the impact parameter bi and 
charged multiplicity N c h (\rj\ < 1), respectively. One sees 
in Fig. 2 (a) that the theoretical results can reproduce 
the STAR data for 0-10% most central collisions. How- 
ever the theoretical results are all higher than the cor- 
responding STAR data for the 10-20%, 20-30%, 30-40%, 
and 40-50% most central collisions and those theoretical 
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FIG. 2: (Color online) The FB multiplicity correlation 
strength b in Au+Au collisions at ^ snn~200 GeV for cen- 
trality determination of (a) impact parameter bi and (b) 
charged particle multiplicity N c h (\r)\ < 1). The open sym- 
bols with line and solid symbols are the PACIAE results and 
STAR data @], respectively. 



results are closed to each other. That is because the FB 
multiplicity correlation in Au+Au collisions is mainly 
the statistical correlation steaming from the multiplic- 
ity fluctuation 14]. And the multiplicity fluctuation in 
those centralities defined by impact parameter are simi- 
lar to each other. 

Figure 2 (b) gives the FB correlation strength calcu- 
lated with the charged particle multiplicity (|ry| < 1) 
centrality determination together with the STAR data. 
We see in this panel that the STAR data are well re- 
produced by the PACIAE model for the 0-10%, 10-20%, 
and 20-30% central collisions. However, the PACIAE re- 
sults are higher than the corresponding STAR data for 
30-40%, and 40-50% central collisions. The correlation 
strength approaches an exponential function of A?/ ob- 
served by STAR in 40-50% central collisions can not be 
reproduced especially. 

In Fig. 3 we compare the FB multiplicity correla- 
tion strength as a function of A77 calculated by the 
wounded nucleon model [ijj], GMC code with a "toy" 
wounded- nuclcon model, HSD transport approach [13j . 
and PACIAE model with the STAR data for the 0-10% 
(panel(a)) and 40-50% central collisions (panel (b)). We 
see in Fig. 3 (a) that all of the four models can nearly 
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FIG. 4: (Color online) The same as Fig. 2 (b) but the PA- 
CIAE results are calculated with STAR's convention of differ- 
ent centrality determination for different measured Ar/ points 
(see text for the details). 



FIG. 3: (Color online) The FB multiplicity correlation 
strength b calculated by different models are compared with 
the STAR dataQ for (a) 0-10% and (b) 40-50% central 
Au+Au collisions at % /sjvjv=200 GeV. 



In this kind of theoretical calculations, each definite cen- 
trality curve in Fig. 4 is composed of b calculated at dif- 
ferent Arj with different centrality determination, such 
kind of curve is not reasonable in theoretical physics. 



reproduce the STAR data for 0-10% central collisions 
and the PACIAE model is somewhat better than oth- 
ers. However all of the model results are higher than the 
STAR data for the 40-50% central collisions, therefore 
can not reproduce b as an exponential function of Arj 
observed by STAR as shown in Fig. 3 (b). 

We have noted the mentions in [6[ that "To avoid a 
bias in the FB correlation measurements, care was taken 
to use different pseudorapidity selection for the centrality 
determination which is also based on multiplicity. There- 
fore, the centrality determination for the FB correlation 
strength for Ar;=0.2, 0.4, and 0.6 is based on the mul- 
tiplicity in 0.5 < \t)\ < 1.0, while for A?j=1.2, 1.4, 1.6, 
and 1.8, the centrality is obtained from \r/\ < 0.5. For 
A7y=0.8 and 1.0, the sum of multiplicities from \rj\ < 0.3 
and 0.8 < \rj\ < 1.0 is used for the centrality deter- 
mination." So we follow STAR's centrality determina- 
tion convention to repeat all of the PACIAE calcula- 
tions and draw Fig. 4 with these results. Comparing 
Fig. 4 with Fig. 2 (b) we see this complicated central- 
ity determination does not improve but even worsens 
the agreement between experiment and theory. That 
means STAR's centrality determination convention may 
be needed for the experimental measurement of FB cor- 
relation strength but not for the theoretical calculation. 
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FIG. 5: (Color online) The PACIAE calculated FB multiplic- 
ity correlation strength b in Pb+Pb collisions at y / 5jvjv=5500 
GeV with centrality determination of charged particle multi- 
plicity N ch (\rj\ < 1). 

Figure 5 gives the PACIAE model prediction for the 
charged particle FB multiplicity correlation strength in 
Pb+Pb collisions at v /saT/v=5500 GeV. We can see that 
the FB multiplicity correlation strength in Pb+Pb col- 
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lisions is smaller than in Au+Au collisions. One reason 
may be that the multiplicity in Pb+Pb collisions at LHC 
energy is much larger than in Au+Au collisions at RHIC 
energy (the N c h/drj at mid-rapidity, \rj\ <0.5, is about 
600 in 0-10% most central Au+Au collisions at RHIC en- 
ergy but it is around 1200 in 0-10% most central Pb+Pb 
collisions at LHC energy [24|). This observation is sim- 
ilar to the report that the elliptic flow parameter V2 in 
Pb+Pb collisions at LHC energy is significantly smaller 
than in Au+Au collisions at RHIC energy [25[. That 
is attributed to the fact that the hard process is more 
influential at the LHC energy than RHIC energy in 25[ . 
Whether the competition between the hard and soft pro- 
cesses is also the reason of the FB multiplicity correlation 
decreasing from RHIC energy to LHC energy is beyond 
this paper scope, and it would be studied later. 



IV. CONCLUSION 

We have used a parton and hadron cascade model, 
PACIAE, to study the centrality dependence of charged 
particle FB multiplicity correlation strength in 0-10%, 
10-20%, 20-30%, 30-40%, and 40-50% central Au+Au 
collisions at ^7^=200 GeV. For the 0-10%, 10-20%, 
and 20-30% central collisions, the STAR data are well 
reproduced. The STAR observed characters of (1) & as a 
function of A77 is approximately flat for central collisions 
and (2) b decreases with decreasing centrality are repro- 
duced as well. However the PACIAE results are higher 
than the STAR data for the 30-40% and 40-50% central 
collisions and can not obtain b as an exponential function 



of At] for the 40-50% central collisions, especially. 

It turned out that the PACIAE model is somewhat 
better than the wounded nucleon model, the CMC code 
with a "toy" wounded-nucleon model, or the HSD trans- 
port model in comparing with the STAR correlation 
data. However all the models can not reproduce b as an 
exponential function of Ar] in the 40-50% central colli- 
sions observed by STAR. That should be studied further. 

The PACIAE calculations are repeated using the 
STAR's centrality determination convention mentioned 
above. The results not improve but even worsens the 
agreement between theory and experiment. This means 
STAR's centrality determination convention may be 
needed for the experimental measurement but not for 
the theoretical calculations. Because in such theoretical 
calculations, each definite centrality curve is composed 
of b calculated at different A77 with different centrality 
determination, this kind of curve is not reasonable in 
theoretical physics. 

A prediction for the charged particle FB multiplic- 
ity correlation strength in 0-10% Pb+Pb collisions at 
Y / s]v r /v=5500 GeV is also given. The charged particle 
FB multiplicity correlation strength in Pb+Pb collisions 
at LHC energy is much smaller than in Au+Au collisions 
at RHIC energy. The further study is out of present pa- 
per scope and has to be investigated in another paper. 
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